Ordering and Melting in Colloidal Molecular Crystal Mixtures 
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We show in simulations that a rich variety of novel orderings such as pinwheel and star states can 
be realized for colloidal molecular crystal mixtures at rational ratios of the number of colloids to the 
number of minima from an underlying periodic substrate. These states can have multi-step melting 
transitions and also show coexistence in which one species disorders while the other species remains 
orientationally ordered. For other mixtures, only partially ordered or frustrated states form. 
PACS numbers: 82.70.Dd 



Recently a new type of colloidal state termed colloidal 
molecular crystals (CMCs) for colloids interacting with 
two-dimensional periodic substrates has been proposed in 
simulations [1,2], observed experimentally [3], and stud- 
ied theoretically [4,5]. CMCs occur when the number of 
colloids is an integer multiple n of the number of sub- 
strate minima, and the n colloids residing in each mini- 
mum form a composite n-mer object such as a dimcr or 
trimer. Neighboring n-mers can have an orientational or- 
dering relative to one another in addition to the spatially 
periodic ordering imposed by the underlying substrate. 
The CMCs were named in analogy with molecular crys- 
tals, in which the molecules have orientational as well as 
crystalline ordering. 

CMCs exhibit interesting two-step melting transitions. 
The first stage of melting occurs when the n-mers begin 
to rotate and lose orientational ordering, while in the 
second stage, the n-mers dissolve and individual colloids 
begin to diffuse throughout the sample [1,3]. The ori- 
entational melting was shown to fall into an Ising-like 
transition class [4,5]. A reentrant disordering can also 
occur in which an orientationally ordered n-mer state 
disorders when the substrate strength is increased [2,3]. 
Recent theoretical work has shown that the orientational 
ordering arises due to the anisotropic multipole interac- 
tions between the n-mers [4]. The reentrant disordering 
occurs when the n-mer is compressed by the increasing 
substrate force, decreasing both the distance between the 
colloids and the multipole moment, until the multipole 
interaction energy drops below the thermal energy. 

When the substrate is varied, new states emerge. Stud- 
ies of colloids interacting with arrays of small traps em- 
bedded in a smooth background have shown coexistence 
between an interstitial colloidal liquid and a pinned or- 
dered solid [6]. In the case of colloids driven over two- 
dimensional (2D) substrates, interesting locking effects 
occur when the colloids move preferentially along certain 
symmetry directions of the underlying lattice [7-10] . Col- 
loidal crystallization on periodic substrates may prove to 
be a useful method to create novel colloidal structures for 
photonic band gap device or filter applications [11]. Ad- 
ditionally, colloids are an ideal model system for study- 
ing collective particle states as well as topological defects 
generated on periodic substrates, which is relevant to a 



variety of other systems including atoms and molecules 
adsorbed on surfaces [12,13], superconducting vortices in- 
teracting with periodic pinning sites [14], or vortices in 
Bose-Einstein condensates interacting with periodic op- 
tical traps [15]. 

An open question is what kind of ordering and melt- 
ing transitions occur for CMCs that are not pure, but 
are composed of a mixture of n-mers and m-mers, where 
ni = n + 1. In this work we show that when the num- 
ber of colloids Nc is a rational, noninteger multiple p of 
the number of substrate minima N^, a variety of novel 
crystalline and partially crystalline colloidal states form 
which cannot occur at integer fillings. These states in- 
clude superlattices of n- and m-mer mixtures for certain 
values of p. As one example, we demonstrate pinwheel 
ordering similar to that observed for molecular dimers ad- 
sorbed on 2D atomic substrates. Other orderings, such 
as star shapes in which dimers point to a monomer, are 
also possible. These mixed states show a rich multi-step 
melting behavior in which the n-mers disorder one species 
at a time, followed by the dissolution of all of the n-mer 
states. The partially crystalline fillings may offer a con- 
venient realization of an orientational glass [16]. 

We consider a 2D system of Nc colloids with periodic 
boundary conditions in the x and y directions, and em- 
ploy Langevin dynamics as used in previous CMC simu- 
lations [1] . The equation of motion for a colloid i is 



dt 



(1) 



Here ff^ = — X^j^i ^^^(''y ) the interaction force from 
the other colloids, where we assume a screened Coulomb 



interaction, V{rij) = {Q/\ri 
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Q = 1 is the charge of the particles, 1/k is the screen- 
ing length, and rj^^) is the position of particle The 
substrate force comes from a triangular substrate such 
as that used in recent experiments [3], and is given 
by = X^Li ^'5in(27r6j/ao)[cos(6'j)x - sin(6'i)y], where 
bi = 2:cos(6'i) - ysin{9i) + ao/2, 9i = 7r/6, 02 = 7r/2 
and 03 ~ 57r/5. Here A is the strength of the substrate. 
We measure length in units of the fixed substrate lattice 
constant oq, and take k — 2/aQ. We have also considered 
the case of square substrates. The thermal force is a 
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FIG. 1. The colloid configurations (black dots) at T = 

for a triangular 2D periodic substrate with A = 2.5. (a) 
A filling of Nc/Ns — 3/2, where a disordered or frustrated 
configuration occurs, (b) A filling of Nc/Ns = 7/4, where a 
star ordering occurs. 

randomly fluctuating force from random kicks. To ob- 
tain initial configurations, we start the system at a tem- 
perature where all the colloids are diffusing rapidly, and 
gradually cool to T = 0. In this model, wc do not take 
into account hydrodynamic effects or possible long-range 
attractions between colloids. We focus on the regime 
where A is sufficiently large that the colloids arc eas- 
ily trapped. In the optical trap experiments, A can be 
readily tuned by varying the laser strength. We have 
checked our resiilts for different system sizes and screen- 
ing lengths and find no qualitative changes. To charac- 
terize the phases, we measure the finite time particle dis- 
placements Sr = (1/Nc) X^iii k«(^)^ 1 1 where t' < t, 
as well as the fraction of six- fold coordinated colloids Pq. 

We first consider mixtures of monomer and dimer col- 
loid states, 1 < Nc/Ns < 2. Previous work has shown 
that a herringbone structure is stabilized at Nc/Ng = 2 
(see Fig. 2(b) of Ref. [1]). At NJN^ = 3/2, illustrated in 
Fig. 1 (a) , we find a disordered arrangement of dimers and 
monomers. The system is frustrated since it is not pos- 
sible to tile the triangular lattice with an alternating ar- 
rangement of monomers and dimers. Thus the monomer 
and dimer locations are disordered, and in addition, the 
dimers have only local orientational ordering. For Nc/Ns 
close to and including 3/2, the freezing temperature for 
dimer rotations is very low due to the lack of orienta- 
tional order. In contrast, at Nc/Ns = 7/4 we find an 
ordered state which we refer to as star ordering, shown 
in Fig. 1(b). Every other minimum in every other row 
captures a monomer while the remaining sites are occu- 
pied by dimers. All the dimers align such that they point 
toward the monomers, producing star like structures. For 
Nc/Ns = 7/4 the orientational disordering temperature, 
where the dimer orientational ordering is lost but the 
colloids remain localized, closely coincides with the ori- 
entational disordering temperature Tc of the herringbone 
state for Nc/Ns = 2. 
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FIG. 2. (a) Colloid configuration (black dots) for filUng 
of Nc/Ns = 7/4 at A = 1.5, where a pinwheel state occurs. 
Two example pinwheels are highlighted with white circles and 
the rotational direction is marked with an arrow, (b) Col- 
loid configuration for an orientationally frustrated filling of 
Nc/Ns = 5/3 at A = 2.5. Local pinwheel ordering is indicated 
by the arrows. An example distorted pinwheel is highlighted 
with white circles. 



The star state shown in Fig. 1(b) only occurs for suf- 
ficiently large A > 2.0. When the substrate is weaker, 
for 0.5 < A < 2.0, a pinwheel ordering occurs instead, 
in which the dimers surrounding the monomer have an 
additional tilt as illustrated in Fig. 2(a). There are two 
possible directions for this tilt and Fig. 2(a) shows the 
counterclockwise case; for different initial random seeds 
there is a 50% probability of producing a clockwise ar- 
rangement instead. Pinwheel states are predicted for the 
adsorption of dimer or linear molecules on triangular sub- 
strates at densities such that a portion of the dimers 
stand perpendicular to the surface, surrounded by the re- 
maining dimers which lie parallel to the surface [12,13]. 
They can also occur for adsorbed atom-molecule mix- 
tures such as CO(i_a:)Arx [17]. Pinwheel states with a 
four-molecule basis are expected for a filling of 7/4 [12], 
and pinwheel states with a larger basis are also possi- 
ble for higher fillings such as 13/7 [13]. For A < 0.5 
at Nc/Ns = 7/4, the colloid- colloid interactions begin to 
dominate over the substrate interaction and the colloids 
form an incommensurate modulated triangular lattice. 

By frustrating the system, we can create an orienta- 
tionally disordered state. For example, in Fig. 2(b) we 
show a partially ordered state that occurs at a filling of 
Nc/Ns = 5/3. Here, unlike the case in Fig. 1(a), there 
is positional ordering, since every third site in every row 
contains a monomer. The dimers, however, do not form 
an orientationally ordered arrangement, although locally 
the dimers have pinwheel or distorted pinwheel struc- 
tures. Both chiralities of the pinwheels appear, as in- 
dicated in Fig. 2(b), but there is no long range dimer 
orientational ordering. At fillings Nc/Ng — 5/4 and 
Nc/Ns = 4/3, the dimers form a triangular lattice; how- 
ever, since the spacing between neighboring dimers is 
larger than ao in these configurations, the 
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FIG. 3. The colloidal positions (black dots) and trajecto- 
ries (lines) over a fixed period of time for A = 3.0 and a 

filling of Nc/Ns = 9/4 at different temperatures, (a) Phase 
I, trimer-dimer solid, at T/T,n = 0. (b) Phase II, trimor 
liquid, dimer solid phase, at T/Tm = 0.56. (c) Phase III, 
trimer-dimer liquid, at T/Tm = 2.4. (d) Modulated liquid 
phase, at T/Tm. = 5.0. 

quadrupole-quadrupole interaction strength is too small 
to impose orientational order and the dimers are rota- 
tionally disordered even at very low temperatures. 

For fillings 2 < Nc/Ns < 3, we obtain mixtures 
of dimers and trimers. At a filling of Nc/Ng = 5/2 
the dimer-trimer arrangement is positionally disordered 
since, similarly to the monomer-dimer case illustrated 
in Fig. 1(a), it is not possible to tile the underlying 
lattice alternately with dimers and trinicrs. At a fill- 
ing of Nc/Ns = 9/4 we find the ordered state shown 
in Fig. 3(a). Here the trimers form a triangular lattice 
where every other substrate minimum in every other row 
captures three colloids. In addition, the trimers are ori- 
entationally ordered, as are the dimers surrounding each 
trimer. 

Since the Nc/Ng = 9/4 dimer-trimer state has two dif- 
ferent types of orientational order, we consider its melting 
properties to see whether the orientational ordering can 
be destroyed in more than one step. We find four phases 
as a function of temperature, illustrated in Fig. 3. We 
report temperature in units of T„j, the temperature at 
which the colloidal lattice at this density melts in the 
absence of the substrate. Phase I, the frozen orientation- 
ally ordered state described above, appears at low T, as 
shown in Fig. 3(a). As T increases, the trimers undergo 



FIG. 4. (a) The average displacement Sr versus temper- 
ature for the system in Fig. 3. (b) Phase diagram for 
Nc/Ns = 9/4 at different values of substrate strength A, 
showing the disordering lines between the different phases. 
Triangles: Tc^\ transition from phase I to phase II; squares: 
Tc^\ transition from phase II to phase III; circles: dissolution 
transition from phase III to a modulated liquid (ML). 

a rotational disordering transition and begin to rotate at 
a temperature Tc^^; however, the dimers remain frozen. 
In Fig. 3(b) we plot the colloid trajectories (black lines) 
over a fixed time period for this phase, termed phase 
II, which shows that the trimers are undergoing rotation 
while the dimers are stationary and ordered. As T is fur- 
ther increased, the system enters phase III at Tc^^ when a 
second orientational disordering transition occurs. Here 
the dimers lose their orientational order and begin to ro- 
tate, as illustrated in Fig. 3(c), but the colloids remain 
trapped within the substrate minima and there is no dif- 
fusion of colloids throughout the system. The trajectories 
of the rotating trimers are widened in Fig. 3(c) compared 
to Fig. 3(b) due to the higher temperatures. If the tem- 
perature is increased sufficiently, the dimers and trimers 
dissolve and diffusion occurs throughout the whole sys- 
tem in a modulated liquid, as shown in Fig. 3(d). Both 
the dimers and the trimers appear to dissolve at the same 
temperature, so we do not find an interstitial liquid of 
colloids moving throughout the sample coexisting with a 
pinned solid, such as seen in systems where the substrate 
has flat regions between traps [6,18]. 

The transitions between the different phases can be ob- 
served by monitoring the average particle displacements 
Sr. In Fig. 4(a) we plot 6r vs T/T„i for the system 
from Fig. 3. Here, Sr = in phase I. The transition 
to the trimer liquid, dimer solid phase II appears as a 
step in Sr at Tc^\ while a second step occurs at tP 
at the transition to the dimer-trimer liquid phase HI. In 
the modulated liquid phase 5r increases rapidly. Note 
that although there are sharp jumps in Sr, these are 
not signatures of a first order transition but instead arise 
due to our definition of the measure. We have used Sr 
to identify the transitions between the phases in a se- 
ries of simulations for different substrate strengths A at 
Nc/Ns = 9/4. We consider only A > 1, since at this 
filling for A < 1 new incommensurate modulated tri- 
angular ordering states appear as the colloid-colloid in- 
teractions become dominant. As indicated in Fig. 4(b), 



the transition temperatures ri^-* and ri^' both decrease 
as A increases, with Tc^^' decreasing more rapidly than 
Tc'^\ In contrast, the crossover hne from phase III to 
the modulated liquid increases linearly with A. In sim- 
ulation [1] and experiments [3] on single-species CMCs, 
the transition temperature from the oriented CMC state 
to the orientationally disordered CMC state decreases 
with increasing substrate strength, as also seen here for 
the multi-species CMC. Theoretical work for the single 
species case [4] has shown that this is due to the reduction 
of the effective multipole moments when the colloids in 
each substrate minimum are forced closer together with 
increasing A. In phase I of the multi-species case, the 
long range orientational order of the trimers arises due 
to interactions with the other trimers. Because the aver- 
age trimer-trimer spacing is much larger than the average 
dimer-dimer spacing, the trimer orientational ordering 
energy is smaller than the dimer orientational energy and 
thus the trimers undergo a rotational disordering transi- 
tion at a lower temperature than the dimers. We note 
that similar multi-step melting transitions occur for other 
filling fractions whenever both n- and m-mer species have 
orientational ordering at low temperatures. 

In conclusion, we have shown that in colloidal molecu- 
lar crystals composed of monomer-dimer or dimer-trimcr 
mixtures, a rich variety of novel ordered and partially 
ordered states can occur. The ordered monomer-dimer 
states include star and pinwlieel states, where the lat- 
ter are similar to phases observed for molecular dimers 
adsorbed on periodic substrates. Unlike the molecular 
system, it is possible to obsc;rvc the colloidal versions of 
these phases directly in experiment, and in addition the 
strength of the quadrupole interaction between dimers 
can be tuned by varying the substrate strength. For other 
monomer-dimer fillings, partially ordered or frustrated 
states can occur, including a positionally ordered but 
orientationally disordered state which may offer a con- 
venient realization of an orientational glass. For ordered 
states in which the two species both have orientational or- 
dering, such as dimer-trimer mixtures, a multi-step melt- 
ing occurs where one species rotationally disorders while 
the other species remains orientationally ordered. As the 
temperature increases further, the second species loses 
its orientational order and finally at high temperatures 
the n-mer states dissociate into a modulated liquid. In 
addition to colloids on periodic substrates, similar states 
may be realizable for vortices in superconductors with 
periodic pinning or for vortices in Bose-Einstein conden- 
sates interacting with periodic substrates. 
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